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This paper is devoted to the study of some aspects of the instability of two flavor color supercon- 
ductive quark matter. We find that, beside color condensates, the Goldstone boson related to the 
breaking of U(1)a suffers of a velocity instability. We relate this wrong sign problem, which implies 
the existence of a Goldstone current in the ground state or of gluonic condensation, to the negative 
squared Meissner mass of the 8 th gluon in the g2SC phase. Moreover we investigate the Meissner 
masses of the gluons and the squared velocity of the Goldstone in the multiple plane wave LOFF 
states, arguing that in such phases both the chromo-magnetic instability and the velocity instability 
are most probably removed. We also do not expect Higgs instability in such multiple plane wave 
LOFF, at least when one considers fluctuations with small momenta. The true vacuum of gapless 
two flavor superconductors is thus expected to be a multiple plane wave LOFF state. 

PACS numbers: 



I. INTRODUCTION 



It is by now well accepted that quark matter at high densities and low temperatures finds itself in a color- 
superconductive state [l| . The ground state at asymptotically large quark densities is the CFL (color-flavor-locked) 
phase 0]. Reviews of high density QCD are given in [H, H, [f| [g, 0, H, H, EH UH EH • It is also admitted however that, 
if the cores of some compact stars contain quark matter, such matter should be in a different pre-asymptotic phase. 
The reason is that the expected values of the chemical potential within such cores would be low enough to make 
non negligible the strange quark mass, requiring a different type of pairing than for CFL. The identification of such 
a phase, respecting electric and color neutrality and weak equilibrium, is an active subject of research at this time, 
and some solutions have been proposed. Homogeneous and non homogenous fcrmion condensates have both been 
examined, the latter case seeming at this time to be better resistant to pathologies. The homogeneous choice leads to 
g2SC in the case of two flavo rs UM . corresponding to analogous situations in condensed matter [l4j]- For three flavors 
it leads to gCFL 0, El E3, M, M, HI ■ 

Both g2SC and gCFL are gapless phases since the fermion dispersion laws may have nodes for some value of the 
quasi-quark momenta (in ordinary superconductors the fermion quasi-particles have node- less dispersion laws). As a 
consequence there exist regions in the phase space where fermion pairings become impossible. The pathology with 
these phases is however a chromo-magnetic instability, both for two flavors [2l]|. and for three flavors (22j. Phases 
with inhomogeneous quark condensates 0, [H, HEj2p, [13, HI [U, H EL [lITHi, H, [H, HI (the so called LOFF 
phases, from condensed matter works in the sixties 37j ) appear as more promising at this time. Beside them, gluonic 
phases [38| and ground states with spontaneous current generation [39|, [if| have been proposed. 

One purpose of the present note will be to sharpen the diagnostic for the chromo-magnetic pathology of the gapless 
phases. We shall add a new element to such a diagnostic by finding that in addition to the negative squared Meissner 
mass, which first revealed the pathology, another element, namely a negative squared velocity of the Goldstone related 
to the axial U(l) current [4f[, is present. Such a Goldstone is not one of the Goldstones to be eaten up by the gluons 
after color gauging. We shall discuss the relation between the two elements: Meissner masses and £^(1)^ Goldstone 
velocity. As the search for the correct pre-asymptotic color superconductive phase cannot be said to be concluded at 
this stage, a full knowledge of the possible pathologies will presumably be useful before concluding such a search. 

We shall write down the low energy lagrangian of 0, the Goldstone boson related to the spontaneous breaking of 
U(1)a- In the gapless regime the Goldstone particle </> is found to exhibit a negative squared velocity. This negative 
velocity will be seen to be directly related to the negative squared Meissner mass of the gluon 8 th . 

The relation with the Meissner chromo-magnetic instability leads to the hope that such Goldstone pathology does 
not apply to the crystalline phases, which arc at this stage the favored candidates for QCD at pre-asymptotic densities, 
beside gluonic phases. We shall examine in the main part of this paper the nature of the true ground state of two 
flavor superconductivity. We shall first verify the equivalence of the one-plane-wave LOFF state with the state with 
Goldstone currents, the cur-g2SC state. Wc shall express the ground state in terms of the two flavor LOFF phase for 
one-plane-wave. We shall then discuss the effective action for the fluctuations of the phase of the condensates in the 
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presence of Goldstone currents, showing that there is no sign problem for the squared velocity. Such a result is shown 
to hold also for multiple-plane-wave LOFF. Contrary to one-plane- wave, such states may be valid candidates for the 
ground state. We discuss the Meissncr masses for multiple-plane- waves, arguing that they are of positive squares. 
Finally we examine the Higgs instability for multiple-plane-wave LOFF, arguing in favor of no Higgs instability at 
least for small momenta. Our conclusion is in favor of a multiple-plane-wave LOFF for the ground state of gapless 
two flavor superconductors. A comparison with the phases of condensed gluons, following the treatment depicted in 
Ref. [13] for the case of the single plane wave LOFF phase, is at this stage of vital importance. 

The plan of the paper is as follows: in Sec. II we introduce the model. In Sec. Ill we derive the parameters of the 
small momenta action for the Goldstone field (j>: the squared decay constant and the squared velocity. In Sec. IV 
we study the true ground state of the two flavor color superconductor. We first report on the equivalence of the 
one-plane-wave LOFF state with the Goldstone current. Then we discuss the fluctuations in the cur-g2SC phase. 
We write down the effective action for the Goldstone boson in the LOFF phase and calculate the Meissner masses in 
the multiple-plane- wave phases. Finally we briefly discuss the Higgs instability for such phases. In the Appendix we 
derive the screening Debye and Meissner masses of a fictitious gauge boson. The calculation is then matched to the 
low energy properties of the </>. 



II. THE MODEL 



In this paper we consider two flavor superconductive quark matter whose action is given by 



d A x 



aB 



■4> jP + (L -> R) + £ A 



(1) 



in the above equation a, f3 = 1,2,3 denote color and i,j = 1,2 stem for flavor. The spinor tjj is a left-handed Weyl 
spinor. The chemical potential matrix fi"f is defined as 



aB 



(2) 



In Eq. @ /ib is one third of the baryon chemical potential; /i e is the electron chemical potential, coupled to quarks 
via the electric charge matrix Qij = Diag[2/3, —1/3,-1/3]; finally /Z3 and /xs are color chemical potentials related to 
the conserved charges Q3 = (ip^T^ip) and Qg = \ip Tstp)- It has been shown in [43| that in the 2SC and the g2SC 
phases of QCD, to which we are interested in this paper, the introduction of /i.3, /^g is enough in order to properly 
achieve color neutrality. 

Condensation in the quark-quark channel is described by the condensation lagrangian £a which is given by 



A rp 



Ct/jjae 



aB3 t 



h.c- (L -> R) 



(3) 



it can be obtained in the mean field approximation from a local four-fermion interaction. In writing Eq. 
assuming that in the ground state 



we are 



(4) 



where the superscripts L, R denote left-handed and right-handed quarks respectively. In the ground state characterized 
by the above quark condensate one has ^3 = 0: as a matter of fact, the vacuum expectation value (j4|) is neutral with 
respect to color transformations generated by T3; therefore Q3 = identically, and the color chemical potential 
associated to T3 vanishes. As for /ig, it has been shown that it is non-zero but in any case it is negligible (that is 
fis "C /i e )- So from now on we will put fig — and assume that the only difference in chemical potential among quarks 
arises from the charge chemical potential (iq = —fi e . 

In this paper we adopt the high density effective theory (HDET) approximation in order to derive the quark 
propagator 0, [44|. This approximation is justified since quarks live at high baryon chemical potential fib and the 
relevant momenta for the dynamics are those near the Fermi momentum (that in the case of massless quarks coincides 
with their chemical potential). In this case, negative energy fields are decoupled and suppressed with respect to the 
positive energy ones. Thus at the leading order in the expansion 1/ ' \i one can describe the system in terms of the 
positive energy fields only. In HDET one decomposes the quark momenta as 



P = (A* + ' 



(5) 
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with v the Fermi velocity of the quarks, \v\ = 1 in the massless case considered here. In the above equation 
M = (m« + Md)/2- The HDET action in momentum space reads 



_ (j^ fdn f di^dlp f f V ■ £ + SfiA 
S ~ it J 8tt J (2tt) 2 X 



(0) 



A V ■ £ - 5nA 

where Sfi = (/id — /i u )/2; x is a left-handed and positive energy velocity dependent Nambu-Gorkov field defined as 

and $ = (ipur, ^ug^dr, ipdg, ipub, ipdb)- In Eq. © we have introduced the matrix A = Diag[— 1, —1, 1, 1, —1, 1]. Finally 
the gap matrix A is 



(8) 



The fermion propagator D(£q,£\\) for left and right handed fields can be easily obtained from the above action once 
we write it in the form 
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n 2 fdn f dindin , , s , 

-j^j -'x + ^V- 



The poles of D define the dispersion law for the quasi-particlcs. They read as 



E±± = ±5{i ± + A 2 , red and green quarks 
E±± — ±(5/x± \£»\ , blue quarks . 

When <5/i > A one can introduce a blocking region (BR) defined as 

{*„ 3' < vV - A 2 } . 



(9) 

(10) 
(11) 

(12) 



III. EFFECTIVE ACTION FOR THE GOLDSTONE FIELD 



In order to derive the effective action for the Goldstone boson </> associated to the breaking of the £^(1)^ symmetry 
in the 2SC and g'2SC phases, we introduce the external field U = cxp(i<p/ f) by means of the following transformation 
over the fermion fields 1451] : 



so that the fermion action becomes 



- r fdn f d^dio vt ( V ■ £ + 6^A Ae 2l ^f 



ir J 8tt J (2tt) 2 X V Aer 2l <f>/f V-i-SfiA 



(13) 



(14) 



The transformation in Eq. (|13p gives rise, at the leading order, to a three body and a four body interaction term 
among quarks and cj>, namely 



j'2£0^ n 2 fdn f dt\\de ^ 



f J 7T J 8tt J (2?r) 



A 
-A 



(15) 



2<t> 2 \ n 2 



dn f d£nd£ Q 



8tt J (2?r) 5 



A 
A 



X + (L -» R) 



(16) 
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Integration over the fermion fields in the functional integral gives rise to the effective lagrangian which, at the second 
order in <fr, consists of two terms, namely a self-energy and a tadpole action: 



2/o 



Tr 



D 



A 

-A 



D 



A 

-A 



(17) 
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Tr 



D 



A 
A 



Evaluation of the traces gives in momentum space 
<j>(-p)<f)(p) p 2 f dn 



iC s . e .(po,p) 



r- 



x / dt\\dt 



and 



4 A 2 [A 2 - (£ + Sp)(£ +p + Sp)+£ ll (£ ll +p- n)} 
[(£ + 5fi) 2 -£\- A 2 ] l(£ +p + Sp) 2 - (£,| + p ■ n) 2 - A 2 ] 



Ctad = ~£s.e.(P0 = 0,p = 0). 



(18) 



+ Sfi -5ft } (19) 



(20) 



The last equation implies that C s , e .(po = 0,p = 0) + Ctad = 0, which is equivalent to the vanishing of the mass of < 
as required by the Goldstone theorem. 



A. Decay constant 



For the computation of the decay constant / it is enough to evaluate £ s , e .(poiP — 0) + Ctad- To this end, in order 
to properly treat the infrared divergences arising from the gapless modes in the fermion propagator, we perform the 
£0 integral at finite temperature, J d£o = 2mT Y] and £q — > iirT(2n + l),po —> iirTm in the integrand. Once the 
summation over Matsubara frequencies is performed we take the limit T — ► 0. The result is 



Ctad + £ s . e .(po,P= 0) 



H-Po) ( f>(Pa) V 2 
P 



X 



(I 



d£ l{ [l-e{5fi-Em))} 



8A 2 



8A 2 



8A 2 



2E{£\\) + Po + i0+ 2E{£ ] \) - p {) - i0+ E{£\{) 



(21) 



In the above relation E(x) = y 1 x 2 + A 2 . In the analytical continuation from imaginary to real boson energy we add 
a small positive imaginary part to pq. As usual, this gives rise to an imaginary part of the polarization tensor for the 
boson, which is related to its decay rate. The imaginary part develops when p$ > 2A, the sum of the rest energies 
of the two quasi-particles. This is equal to what happens in the gapped 2SC phase, when Sfi = 0. The integral over 
quark momentum can be performed analytically for each value of po; at the order j»q/A 2 one gets 



4^ 2 



£tad+ £ S .e.(P0,P = 0) = ^ (PO<£)(PO0)-^2- I 1 - 0(Sp - A) - 



(22) 



Imposing the canonical normalization of the field <fi in Eq. 



fixes the value of /, namely 



f = V 1 _ e{ s ll _ A) 



ysp 2 - a 2 

Sfi 



(23) 



When SfL = we obtain the well known result of the 2SC phase, f^sc = 4/x /7r . In the limit Sp^> A one has 



2_/^ A_ 2 

7T 2 8fl 2 



f ~ _9 !..0 



(24) 
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B. Squared velocity 

Next we turn to the computation of the squared velocity of <f>. To this end it is enough to consider C s . e . (0, p) + Ctad- 
The integral over quark energy is performed as above at finite temperature, keeping the limit T — > after the 
summation on the Matsubara frequencies is done. We find 

C s . e .(0,p)+C tad = (t>{ -^f p) !l |fg J +O °d£ {] [F(l h p-n)-F(£ h 0)] (25) 

with 

F(x, y) = ; 2X ( 1 - 0(5fx - Vx 2 + A 2 )) + 2(x + ^ f l _ 0(5^ - \/(x + y) 2 + A 2 )) . 

(26) 

Once the integral over the quark momentum is performed (using the 6— functions) we expand at 0(p 2 / A 2 ); we find 
£ s . e .(0,p) + £ tad = -^(p^)(p^)^(l-0(^-A) 



{Pi<t>)(pi4>) , (27) 



v 2 



2 



with the squared velocity defined as 



and f 2 given by Eq. For Sfi = we recover the well known result v\ sc = 1/3. With these results at hand we 

can write the effective action of d> at small momenta as 



m = J d A x l - {{d 4>) 2 - v 2 {V^V4>)) . (29) 



Eq. |29|) . with / 2 and w 2 given respectively in Eqs. |23|) and (|28|) . is one of the results of this paper. 

Notice that in the gapless phase we find a negative squared velocity of (p. This instability in the Goldstone sector 
is directly related to the chromo-magnetic instability of the gapless 2SC phase. As a matter of fact, we see that we 
can write the relation 



v 2 oc m 2 M 8 (30) 
where m 2 M 8 is the squared Meissner mass of the 8 th gluon calculated in [2l[ 

mlA .^- m -^L^ . ,3!) 

As far as the gapped phases, the relation between the squared velocity of the Goldstone mode and the squared 
Meissner mass of a gluon (or, more generally, of a spin 1 gauge field) is not new, see [H, [HD, US Hil for 2SC and 
CFL. It is interesting to notice that the v 2 of <fi is negative only in the gapless regime Sfi > A, while the 2SC 
phase presents the chromo-magnetic instability also in the gap ped region A/^/2 < 5fi < A: in this interval one has 
m M 8 > but m 2 M a < 0, with a = 4, ■ • ■ , 7. As shown in [38} the instability in the gapped regime can be cured by 

a gluonic condensate. The method of Ref. [38| however can be applied only in the region Sp, rts A/\/2, where the 
gluon condensates are expected to be small, because it relies on the expansion of the expectation values of some of 
the gluonic fields around zero, see [HI for more details. In the gapless regime Sfj, > A the instability could be cured 
cither by a gluon condensate or by net baryon and/or meson currents, see (39j for a wide discussion. We expect that 
the removal of the chromo-magnetic instability for the 8 th gluon is enough for the removal of the Goldstone mode 
instability. Finally, we notice that v 2 < for the Goldstone mode <p is similar to the instability found in Ref. for 
the 8 th gluon. 
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IV. THE GENUINE GROUND STATE 



In the previous section we have computed the low energy parameters of the Goldstone mode <f> both in the 2SC and 
in the g2SC phase. In particular, we have found a negative squared velocity, a signal of an instability of the ground 
state. Therefore, the g2SC can not be the true vacuum of the model and one has to look for other solutions. In 
this section we introduce an ansatz for the true ground state, which imitates the LOFF phase of condensed matter. 
This topic is not new: both in the two flavor and in the three flavor gapless phases of QCD it has been extensively 
studied (2^1 . However, it has never been related to the negative squared velocity of one (or more) Goldstone bosons. 
We cover this topic in this section. 



A. The one plane wave LOFF state and its equivalence with the Goldstone current 

The starting point is the observation that, because of the wrong sign in the spatial part of the kinetic term in 
Eq. ([29]) . the expectation value for (VEV) may be expected to be different from zero. In order to determine the 
true ground state we decompose the <j> around a mean field part, and study the effective action for the fluctuations 
around this mean field. Moreover, since V<^> is a vector, the VEV breaks spontaneously the rotational symmetry. The 
simplest way to introduce the fluctuations is writing <j> as 



4>(x) = $ x + h(x) , 

where $ is a constant and homogeneous vector and h is the fluctuation field; in this way one has 

= * + Vh ; 



(32) 



(33) 



assuming (V/i) = we are left with (V0) = <&. Then, choosing 3? appropriately, one has to show that the lagrangian 
of the fluctuations h(x) does not suffer of the sign problem in the spatial part. Since *& ^ corresponds to a Goldstone 
current, we call the ground state cur-g2SC phase. 

Introducing Eq. (f3"2"|) into Eq. (jT5J) one notices that the action resembles that for the two flavor LOFF phase in the 
one plane wave (1PW) structure, 



S 



d 4 a 



dn 

8^ 



A exp 



iV ■ d + S^iA 
-2i(* • x 4 



A exp 



h) 



2 if* 



f 



iV-d 



f 

SfiA 



(34) 



/ 



with <&// = q playing the role of the wave vector and the field h being the phonon [28|] as well as the Goldstone boson 
related to the U(1)a- Since the fermion dispersion laws in 1PW and in cur-g2SC are the same, the thermodynamic 
behavior of the two phases is indistinguishable and in what follows we can refer both as the 1PW and as the cur-g2SC. 
However we stress that the cur-g2SC phase is built with zero momentum Cooper pairs, and the wave vector q in this 
context is related to the non-vanishing Goldstone current in the ground state. On the other hand, the 1PW ground 
state is built with Cooper pairs with a total momentum equal to 2q; the Goldstone current in this case is vanishing, 
as we prove in the next section. 

The 1PW state can be analyzed exactly, and the effective action for the fluctuation h can be determined without 
approximations: it is enough to shift the quark momenta by an amount ±q, where the upper and lower signs stem 
respectively for u and d quarks. In this way the x dependence in the gap term is ruled out in place of a shift 
6/j, — > Sfj, — q ■ n in the quark chemical potentials; but in this case the expressions involved in the calculations of 
the loop integrals are more complicated than the homogeneous ones. Therefore we prefer to treat this problem by a 
Ginzburg-Landau (GL) expansion of the fermion propagator in A/5/i. In expanding a la Ginzburg-Landau we pay 
the price of loosing exact expressions, obtaining equations valid only to a fixed order in A/Sfi, but the formula are 
easy to handle and to generalize to the case of more complicated crystal structures (see below). 

In the GL approximation one formally writes the quark propagator, which can be read from Eq. ([6]), as 



A 
A 



D , 



(35) 



where Dq is the pro pag ator at A = 0. The zero temperature thermodynamic potential for the quarks in this 
approximation reads [29( 



n = n, 







2 A ' + I A 



(36) 
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where the coefficients are given by 



Q 



4^ 2 



Sfi 



8(i 



1 A 2 

2 1 ° g 4(|q| 2 -<5 At 2 ; 



(37) 







8(i 2 



ir 2 \q\ 2 - 5(i 2 



8tt 2 (|q| 2 - 5(1 



2\3 



(38) 



and Q Q is the free gas contribution. The physical value of q is obtained as usual by minimization of the thermodynamic 
potential. At the leading order in A/8(i one has 



d\q~\ 



q\=Q 



0^1-^log 



Q — Sfi 







(39) 



which gives the result Q ~ 1.2<5/i well known in the LOFF literature. Since 3> 2 = f 2 q 2 , once we know f 2 we are able 
to evaluate the Goldstone current in the cur-g2SC state. 

B. Effective action of the fluctuation in the cur-g2SC state 

Next we turn to the effective action for the fluctuation field h. Evaluating the traces in Eqs. (fTT)) and (fTS)) at the 
leading order in A/8(i one is left with the expression 



£s.e.(p) + £ 



tad 



2A' 



H-p) U(p) - J(0)] h( P ) 



(40) 



where p = (po,p) and the loop integral J(p) is defined as 

d 4 i 



" dn 



4tt J (2tt) 4 



1 



(V -£-S(i + q- n)[V ■ {£ + p) - 8(i + q ■ n] 
1 

(V ■ £ - 8(i + q ■ n)[V ■ (£ +p) - 8(1 + q ■ n] 



5(i 



-5(i . 



(41) 



From now on the calculation is similar to the one presented in great detail in [47j for the displacement fields in the 
three flavor LOFF phase of QCD; therefore in this paper we simply show the main steps of the calculation, referring 
the interested reader to Ref. [47J for further details. For small external momenta p one has 



J(p) ~ J(0) 



dn 



d 4 £ 



V-pV -p 



4tt J (2tt) 4 (V ■ £ - 5(i + q ■ n) 2 {V ■ £ - 8(1 + q 



+ 5(i — ► — 5(i 



(42) 



The computation of the loop integral is done in the usual way by Wick rotating to imaginary energies £q — ► i£i, since 
the integral is convergent one can send the ultraviolet cutoff on £u to infinity, and perform the integral over £\\ by 
residues, followed by integration over £4. This is the same procedure used for the calculation of the coefficients f3, 7 
in the GL effective potential [2^]. We find 



J(P) - J(0) = -i£fle I V ' PV -' +8(1^ -5(i 

4tt 2 J 4?r (8(i-q-n + 1O+) 2 



(43) 



From Eqs. (|4"0"1) and (|4^|) one can easily read the low energy parameters of the effective lagrangian for the fluctuation 
field in the cur-g2SC phase, namely 



47r (S(i — q ■ n + i0 + )'- 



+ (8(i — > —8(1) 



2(i 2 



7T 2 Q 2 - 5(1 2 



(44) 



2 ,,2 



A 2 (i 



■5Re 



dn 



4tt (5(i q n + i0 + )' 



+ (5(i -> -8(j) = 



(45) 
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A 2 /i 2 fdn nl , t 2n 2 A 2 . . 

«. = -JW^J 4^-g.n + l 0+) 2 + = /^Q 2 -^ ■ ( 46 ) 

Finally, in configuration space, the lagrangian of the fluctuation reads 

C[h] = \ ((d h) 2 - v ■ (Vh) v ■ (Vh)) . (47) 
By means of Eq. (|44[) we determine the value of the Goldstone current which minimizes the effective potential, 

|(V^)| 2 = |*| 2 = ^A 2 ^^ « 0.66/, 2 A 2 , (48) 

where we have used the relation Q ~ 1.25 /j,. Moreover Eqs. (|45|) and ([46]) show that the squared velocity of h is positive 
along the direction of and is zero in the plane orthogonal to 3>. As can be easily shown by direct calculation, the 
zero value of the orthogonal velocity is due to the proportionality relation between the integral in Eq. (|45[) and the 
derivative dVt/d\q\ evaluated at the minimum, see Eq. (|39|) . Therefore, the fluctuation h does not suffer the sign 
problem of the squared velocity, as anticipated. 



C. Effective action of the Goldstone boson in the LOFF phase 

In the previous section we have shown that assuming the existence of a Goldstone current in the ground state of 
the gapless 2SC quark matter, and expanding the Goldstone held around the VEV as in Eq. (|32p . the lagrangian 
of the fluctuation does not suffer the sign problem of the squared velocity. The cur-g2SC phase resembles the one 
plane wave LOFF state in the sense that the breaking of the translational and rotational symmetries due to <fr 7^ 
is not distinguishable from the breaking due to a net momentum of the Cooper pair, as can be seen by Eq. (|34|l . 
Therefore, it is obvious that the effective lagrangian of the Goldstone of U(1)a in the 1PW LOFF phase is equal to 
the lagrangian of the fluctuation found in the previous section. In the 1PW phase the Goldstone does not suffer of 
the wrong sign problem, and the free energy of the ground state with Goldstone current is equal to the free energy of 
the 1PW phase. 

In the 1PW phase the gap parameter has an explicit spatial dependence of the form 

A(r) 1PW = A e 2i « r , (49) 

with q being the wave vector (we have dropped for simplicity the color, flavor and Dirac indices), equal to one half 
of the total momentum of the Cooper pairs. Because of the anisotropy, only a fraction of the Fermi surfaces of the 
quarks are available for the pairing: this results in A (and consequently, f2) smaller than the one of the BCS phase. 
However the free energy of one plane wave LOFF state can be lowered by summing up P plane waves (PPW) [29j, 

p 

a=l 

the resulting phase is known as a crystalline superconductor, as the behavior of the gap in the configuration space 
resembles that of a crystal lattice. 

Since the free energy of the LOFF phase with order parameter given by Eq. (|50p is lower than the single plane 
wave one, it follows that the 1PW can not be the ground state of a crystalline superconductor. Stated in other words, 
since the cur-g2SC and the 1PW phases are not distinguishable, the free energy of the cur-g2SC phase is higher than 
the free energy of a PPW crystalline state. Therefore it seems that the cur-g2SC can not be the ground state of two 
flavor quark matter, as it can be easily replaced by crystalline phases. This situation is quite different from the three 
flavor case. As a matter of fact, it has been shown in [40( that a Goldstone current there exists in the ground state, 
near the onset CFL— >gCFL. The curCFL phase considered in is not likely to be replaced by a multiple plane 
wave LOFF state near the onset since the free energy of the latter, as evaluated in (3f|, is higher than the gCFL one 
and therefore still higher than the energy of the curCFL phase. 

Since in the crystalline LOFF phases there exists the Goldstone related to the breaking of U(1)a, the calculation of 
its effective action in the PPW state becomes of vital importance. Nevertheless, we learn from the 1PW an important 
lesson: the effective lagrangian, at least in this simple case, does not suffer the sign problem. We wish to verify that 
the same property is valid for a multiple plane wave crystalline superconductor, which is a better candidate for the 
ground state. 
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In the PPW phase the inverse quark propagator can not be inverted, so one is forced to make some approximation 
in order to write a propagator. As anticipated in the previous section we employ the GL expansion. In this case, at 
the leading order in A/Sfi, the value of \q\ which minimizes the thermodynamic potential is again given by Eq. (|39|) . 
in which now Q represents the equilibrium value of the total momentum of the pairs [29| (in the cur-g2SC context it 
is related instead to the VEV of the Goldstone current ) . 

The calculation of the effective action can be done following the same steps outlined in the previous section, replacing 
the fluctuation h with the Goldstone <fi. As shown explicitly in [47j], in the PPW state Eq. (|4"3")) has to be replaced 
with 

J(p)-J(0) = -i£*teY l V-pV-p § s (51) 

4-7T 2 47T (Sfi - q a ■ n + t0+) 2 

In the above equation the interference terms that mix different q a do not appear: this is due to momentum conservation 
in the loop integral, which forces the wave vector of the two gap insertions to be equals. From Eq. (|5Tj) one has for 
the squared decay constant 

PA2 (52) 
with Q ~ 1.2<5/i. For what concerns the squared velocity it is convenient to introduce the matrix 



[ n '" 3 . n+ , 2 + (6ft -> -5ft) ; (53) 

! tt 2 ^ J 4tt (5 ft - q a ■ n + iQ+) 2 v ^ ^' y ' 



the elements of the matrix V are easily evaluated: following the strategy depicted in [30J, for each angular integral in 
the sum one rotates the q a along the z— axes by means of the orthogonal matrix R a . In this way one is left with the 
expression 

1 p 

V « = pYl R li R % . (54) 

a=l 

where we have used Eq. (|52p and the fact that at the minimum only the longitudinal integrals are not vanishing, see 
Eqs. (|45|) and (|46|) . With this definition the effective action for the Goldstone reads, in momentum space, 

m = \<t>{-p) (pg - V« PiPj ) <j>(p) . (55) 

The dispersion law E(p) of the Goldstone is the solution of the equation 

E{pf - Vij p iPj = . (56) 

As a consequence, in order to show that the lagrangian does not suffer the sign problem, it is sufficient to show that 
the matrix V is semi-definite (or definite) positive, that is for each p ^ it satisfies the condition Vij PiPj (or 
Vij PiPj > 0). This is an easy task: as a matter of fact, from the very definition of the rotation matrices R a it 
follows that one can write Vij = J2 a QiQj/Pi w ^h ot a unit vector parallel to q°; since the tensor qf q'- is semi-definite 
positive, as it has one eigenvalue 1 and two eigenvalues 0, then Vij is semi-definite or definite positive, depending on 
the particular crystalline structure considered. We show this in three interesting cases. 

In the first case we consider a crystal with three mutually orthogonal wave vectors (3PW). Each of these vectors 
can be chosen along one of the three axes, 

9l = (3(1,0,0), g 2 = Q(0,l,0), 93 = 0(0,0,1) . (57) 

The second structure we consider is the body centered cube crystal (BCC), defined by the following wave vectors: 

Ql = Q(l,0,0), g 2 = Q(0,l,0), q 3 = Q(0,0,l), 

94 = -9i , 95 = -92 , 96 = -93 • (58) 

Finally we consider the case of the face centered cube structure (FCC), whose wave vectors are given by 

9i = ^=(+1, g 2 = ^=(+l,-l,+l), 
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g 3 = -^L(-l,-l,+l), g 4 = ^L(-l,+l,+l), 
95 = -^(+1+, 1,-1), 96 = -^=(+1,-1,-1), 
97 = ^(-1,-1,-1), 98 = ^(-1, +1,-1) 



(59) 



The rotation matrices R a are trivial for both the case of the 3PW and the BCC crystals. As for the FCC they can 
be found in the appendix B of [3(j. For all of the crystals we find Vy = <5y/3, which implies that the dispersion law 
of the Goldstone is 



E\p) 



1 



(60) 



and the squared velocity does not suffer the wrong sign problem. 

From the results of [29|, [32| we can infer that the BCC and the FCC may exist in a large window of the ratio S/j,/Aq, 
both on the left and on the right of the Clogstone limit Ao/y/2- Since our result shows that in these crystalline phases 
the Goldstone does not suffer the instability towards the formation of a current, we argue that such phases could be 
the genuine ground state of two flavor quark matter at high density. A more accurate study is needed at this point: 
indeed, it would be interesting to build electrical and color neutral crystalline phases in two flavor quark matter, and 
to compare their free energy with the gluonic phases and/or the gapped 2SC. 

D. Meissner masses of the gluons in the PPW LOFF phase 

In this subsection we compute the Meissner masses of the gluons in the multiple plane wave state. The calculation 
is done in the Ginzburg-Landau approximation at the second order in A/<5/i. The Meissner tensor can be defined 
as @ 



0<T 



-n&(po = o,p-o) , 



(61) 



where the polarization tensor in the HDET approach is given by the sum of two contributions, n^(p) ~ S^(p) , 
where S^(p) is a self-energy diagram, 



. 2/i 2 f dn f dl\\dl 



and is a tadpole diagram, 



-Tr 



D(£ + p) 



-V»T a 

o wt: 



D(£) 



-V v T b 







rr\a> 
1 ab 



.2x2 fdn f dt\\dt 0r ^ 



D(£) I 2 "+ v - i 




T a T b 





(M+l||) 2 rpmrp* 



V v Ti 



p„ 



(62) 



(63) 



In the above equations D(i) is the fcrmion propagator, which will be evaluated in the following in the GL approx- 
imation. T a is the ££7(3) color generator in the basis VP = (i^ U r,'<Pug,'4 l dr,'<Pdg,'<Pub,'^db), with the normalization 
Tr[r a T(,] = 6 a b/2. Finally the projector is defined as 



2 



(64) 



The trace in the diagrams have to be evaluated in NG as well as in color-flavor indices. 

Technically, the evaluation of the loop integrals for the Meissner masses is similar to that involved in the calculation 
of the decay constant and the velocity of the Goldstone boson: for all of the gluons one is left with two kinds of 
integrals, namely 



a=l 



El Ull 



dn f d£\\d£o 



4tt ./ (2tt) 2 (V • £ - Sfj,) 3 (V-£-Sfj, + 2q a ■ n) 



+ 5[i — ► — Sfj, 



p r 
-JJeV / 



dn 



167T - X J 4tt (Sfi — q a ■ n + i0+) 2 



+ (5(i — > —Sfi) 



(65) 
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related to a diagram with two gap insertions on a fermion branch and zero insertions on the other branch, and 

related to a diagram with one gap insertions on a fermion branch and one gap insertion on the other branch. Because 
of the proportionality relation between Jy and Jij , one can express the Meissner tensor for each crystal structure only 
in terms of I,. At the leading order in A/Sfj, wc find 



(M^=(Ml) 2 =(A4) 2 = 0, (67) 
(MUf = -i^f-Iij , a = 4,..., 7, (68) 
(M«) 2 = ^(Mii) 2 . (69) 

The above results are in agreement with Eqs. (93), (94) of Ref. [48[, where they have been deduced for the first 
time in the single plane wave LOFF state and have been generalized here to the case of a generic crystalline structure 
(for simplicity here we have neglected the mixing of the 8*^ gluon with the photon. If one keeps into account this 
effect, only the overall factor in Eq. (|69|) changes, as shown in [HI). Looking at Eqs. (|53[) and ([68)) . (|69|) it is easy 
recognized that the following identities hold, 

{Mlf = Qvij , (Mg) 2 = Qv tJ , (70) 

where Vij is the coefficient that multiplies the gradient term in the action of the Goldstone boson U(1)a- Therefore 
the positivity of the matrix Vij reflects to the positivity of the Meissner tensor. Since we have shown that both in 
the case of the BCC and of the FCC crystals Vij is indeed defined positive, a real value for all of the Meissner masses 
follows at the leading order in A/S/i. 

We may ask how the higher order corrections can modify the result expressed in Eqs. (f6"5| . ([59")) . To this end it is 
enough to investigate on the typical ratio A/5/j, in a crystalline LOFF state. Unfortunately, as already stressed in 
the previous sections, the free energy (and the gap parameters) in a generic crystalline structure can be computed 
only in some approximation. In Ref. [32j | a smearing over the cell has been introduced: in this scheme it was found 
that for the BCC the ratio A/Sfi lies in the interval 0.35 - 0.40; for the FCC the interval is 0.27 - 0.30 (see Tables 
I and II of [13]). Since the corrections to Eqs (j68j) . ([69]) are of order A 4 /<5//, we see that the next-to- leading order 
terms are suppressed as A 2 /5^l 2 when compared to the leading order results. At the worst, one has A 2 /8/j? ~ 0.2 
for the BCC and A 2 /5/i 2 ~ 0.1 for the FCC. Therefore we argue that the next-to-leading order terms should not 
change in a dramatic way the main result of Eqs. (f6"5| . ([59")) , namely positive squared Meissner masses for all of the 
gluons in a crystalline LOFF phase. Of course this hint should be supported by the calculation of the next-to-leading 
order correction: unfortunately the computation of such terms in BCC and/or FCC is much more involved than in 
the simple plane wave case. Another strategy would be the calculation using the smearing approach [32| instead of 
the GL expansion of the Meissner tensor. We will come back to this problem in the future. 



E. A brief comment about Higgs stability of the LOFF phases 



Before concluding this section we briefly discuss how the Hi ggs instability could arise in the LOFF phases of more 
than a single plane- wave, following the treatment discussed in [49(. The Higgs (or amplitude) instability is related to 
inhomogeneous variation of the amplitude of gap parameter A. It has been studied in great detail for the homogeneous 
g2SC phase 0, \E^, [5l| , where it has been related to the existence of the Sarma instability. 

In the LOFF phase, one can introduce a variation of the gap parameter in Eq. (|50[) by means of the replacement 
A — > A + H(r) with H(r) real, with Fourier transform given by H(k): 

(71) 

Keeping into account fluctuations in the magnitude of A and neglecting the phase fluctuations (i.e. the Goldstones) 
one can write the grand potential, in the Gaussian approximation, as 

1 f d 3 k d 2 VL 

fi = ^ + 2 / H ^ k) dH^mk) H{k) ' (72) 



Mr) 



PPW 



,i2q a v 



A 



d 6 k 
(2^)3 



-ik-r 



H(k) 
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where fio is the grand potential evaluated at the minimum and the derivative, evaluated at the minimum of the 
free energy, is nothing but the self-energy of the fluctuations (see (Hoj for more details). The momentum integral in 
Eq. (|72|) can be divided grosso modo into an integral over small momenta (with respect to q) and an integral over 
larger momenta. Following [HUANG], in the low momentum region one can replace, in the thermodynamic limit, 

8 2 fl B 2 VL 

(73) 



8H*(k)dH(k) OA 2 ' 

the derivative on the r.h.s. of the above equation is positive because the LOFF state is in a minimum of the grand 
potential. Therefore a LOFF phase can not suffer Higgs instability, at least for small momenta. 



V. CONCLUSIONS 



In this paper we have computed the low energy properties of the Goldstone mode <fi both in the gapped and in the 
gapless 2SC phase of QCD. We stress that such a Goldstone is related to the axial U(l) current and it is not one of 
the would-be Goldstones to be eaten up after gauging of color. In the calculation presented in Sec. IIIII we introduce 
the 4> as an external field. The integration over the fermion fields in the functional integral allows for the computation 
of the squared decay constant and of the squared velocity of <fi. In particular, we find a negative squared velocity in 
the gapless regime 5^ > A. The simple proportionality relation between v 2 and the squared Mcissner mass of the 8 th 
gluon shows that the two instabilities, in the Goldstone and in the gluon sectors, are related and the removal of the 
second is equivalent to the removal of the first one. 

Because of the wrong sign of the squared Goldstone velocity, the gradient of the Goldstone field may take on a 
non-vanishing expectation value in the vacuum, with consequent breaking of the symmetry under rotations. We call 
this phase cur-g2SC phase. The Lagrangian for the fluctuations can be constructed in the cur-g2SC phase and does 
not suffer of the squarcd-velocity sign problem. 

Since the thermodynamics of the cur-g2SC is the same as that of the one-plane-wave LOFF state, its free energy 
is higher than the energy of a multiplc-plane-wave phase. Motivated by this observation we have computed the 
low energy parameters of the effective lagrangian for the U(1)a Goldstone in the multiplc-plane-wave LOFF phase 
(PPW): we find that that such crystalline phases do not present instabilities towards formation of currents. Moreover 
we have calculated the Meissner masses of the gluons in the PPW phases, generalizing the results of [48| : we find 
chromo-magnetic stability. 

The other instability which could appear in in these LOFF phases is the Higgs instability, which has been studied 
for the homogeneous g2SC and related to the Sarma instability. We have argued that the Higgs instability should be 
absent in the PPW LOFF phase, at least for small momenta. A complete study of the amplitude instability requires in 
addition, in the Gaussian approximation, the calculation of the self-energy of H(k) for each value of the 3-momentum. 

A multiple-plane-wave LOFF seems in conclusion to be one of the best candidates for the ground state of a two 
flavor color superconductor, beside the gluonic phases considered in (38|. A comparison of the free energies of the 
neutral LOFF phases in the multiple plane wave state and of the phases with condensed gluons is therefore crucial 
for a deeper understanding of the phase diagram of two flavor QCD (in the single plane wave case this study has been 
performed in [42|). 

Beside the comparison of the PPW state with the gluonic phases, a natural prosecution of this work should be the 
extension to the three flavor superconductive phases of QCD. In this context a lot of work has been done in order to 
remove the chromo-magnetic instability of the gapless CFL phase, see (39[. Moreover, the three flavor superconductive 
crystalline phase does not seem to have such an instability, at least when the Meissner masses are computed in the 
Ginzburg-Landau approximation [52j . It would be interesting to study the Goldstone properties both in the kaon 
condensed and in the three flavor LOFF phases of QCD, in order to see whether the instability found here is present 
also in these cases, and if it is completely removed once the gluon sector is cured. 
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APPENDIX A: AN ATTEMPT TO MATCH LOW AND HIGH ENERGY THEORIES 



In the previous sections we have computed the squared decay constant f 2 and the squared velocity v 2 of the 
Goldstone boson (j>, both in the gapped and in the gapless 2SC phases of QCD. The results are in Eq. (|2"3")l and 
Eq. (|28[) . In particular, we find a decreasing decay constant as the ratio Sfi /A is increased in the gapless phase. 
Moreover, when Sfi > A we find a negative squared velocity of <j). 

For the gapped 2SC and CFL phases it is possible to compute the low energy parameters directly from the screening 
masses of the gluons, or more generally of spin 1 gauge fields, see e.g. [H3, 154L I5H [56^ . We wish to perform such a 
program also in this case, trying to relate the decay constant and the squared velocity of <j) found in the previous 
section to the Debye and Meissner masses of a fictitious gauge boson related to the U(\)a symmetry. In this section 
we compute these masses. 

To begin with we promote the global U(1)a symmetry to a local one: this is done in the usual textbook way by 
the introduction of a fictitious gauge field Wp whose coupling to the fermions is described via the covariant derivative 
Dp = dp + iQaWp- Here Qa is the charge of the fermions under U(1)a, Qa = ±1 respectively for left and right 
handed quarks. The calculation of the static screening masses is similar to that presented in the main text for the 
low energy parameters of the Goldstone mode; therefore here we simply quote the main results. In this case one has 
to evaluate the polarization tensor Up V of the gauge boson Wp. According to Ref [|[ we define 

m 2 D = n o(go = 0, q -> 0) , m\ l S l] = -Hy (q = 0, q -»• 0) . (Al) 
First we consider the Debye mass. After some computation we are left with the expression 

m 2 D = ^-A 2 + M 2 lue , (A2) 



7T 2 J Q (£| + A 2 ) 3 / 

where the first addendum on the r.h.s. is the contribution of the paired red and green quarks, while the second 
addendum is the contribution of the unpaired blue quarks given by the standard many-body result 

M 2 ue = 2x4- (A3) 

In the gapless phase Sfi > A it is interesting to distinguish, in the loop integral involving paired quarks, the contribution 
of the quarks whose loop momentum In lies in the blocking region (|12p from the contribution of the quarks living in 
the pairing region (defined as the complementary of the blocking region): 



+00 




(e 2 + a 2 ) 3 / 2 



(A4) 

where BR and PR stem respectively for Blocking Region and Pairing Region. With these definitions at hand one can 
write 

m% = M% R + M PR + M 2 ue =6x4- (A5) 

7T 

The Debye mass is independent of the value of the ratio Sfi/ A and is equal to the result of the normal phase. This 
can be understood by noticing that both the condensate and the fermion fields have U(1)a charge, so one expects 
screening both in the gapped and in the gapless phase. 
Next we turn to the Meissner mass. We get 

In the gapped 2SC phase one has Sfi < A and the squared Meissner mass is positive and constant. On the other 
hand in the gapless phase Sfi > A the squared Meissner mass is negative and divergent ad the transition point. The 
divergence is related to the divergence of the density of gapless quasi-particle states. 
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The calculations of the Dcbye and Meissner masses which lead to the results in Eqs. (|A5[) and (|A6|) are done using 
the so called high energy theory [54| , that is using a lagrangian written in terms of the quarks degrees of freedom only. 
It would be interesting to make the same calculation in the so called low energy theory, that is by using a lagrangian 
written in terms of the light degrees of freedom, the Goldstone mode and the gapless quarks. The matching of the 
results obtained by the two procedures should allow to compute f 2 and v 2 . 

In the gapped CFL phase of QCD the matching procedure is easy to perform since the light degrees of freedom are 
the Goldstone bosons (an octet plus a singlet), while the high energy theory is defined in terms of the nine gapped 
quarks. Also in the 2SC phase the paired quarks are gapped, so the low energy degrees of freedom are the singlet <f> 
and the unpaired blue quarks. In the gapless 2SC phase, on the other hand, one has also gapless paired fermions in 
the low energy spectrum, so one should consider them in the matching procedure. 

Let us consider first the case of the gapped 2SC phase. On the basis of invariance under U(1)a h is easy to recognize 
that the effective lagrangian for U = exp(icj>/f) can be cast in the form 

f 2 

£ = y (d Ud Q U* - v 2 VUVU*) (A7) 

(with this definition the kinetic term of the field <fi gets the canonical normalization). The gauging of the U(1)a is 
made by the replacement of the usual derivative with the covariant one = + iWp. After the gauging, a mass 
term for the fictitious gauge field is obtained, namely 

SC = (W 2 - v 2 W ■ W) = ^W 2 - ■ W , (A8) 

where we denote by the capital letter Md,m the contribution of the Goldstone to the Debye and the Meissner mass of 
Wp. However one has to add the contribution of the blue unpaired quarks to Eq. (|A7[) in order to properly describe 
the low energy theory. The blue quarks contribute only to the Debye mass, as they are not superconductive and do 
not screen static "magnetic" fields W. As a consequence, the Debye and Meissner masses in the low energy theory 
are 

m 2 D = M 2 D + M 2 lue , m 2 M = M 2 M . (A9) 

On the other hand one can read m 2 D and m 2 M from the high energy calculations, see Eqs. (|A5[) and (|A6[) . Comparison 
allows to compute f 2 and v 2 , 



!lsc = K, v\ sc = \. (A10) 



Now we turn to the gapless 2SC phase and consider the Debye mass. In this case the effective lagrangian for <j> has 
the same form of the gapped phase, Eq. (|A7|) . since it is related only to the symmetries of the ground state and not to 
the absence/presence of gapless excitations in the spectrum. Therefore the contribution of <p to the screening masses 
is described by Eq. (|A8[) and m 2 D is given in the low energy theory by 

m 2 D = M 2 D + M 2 ue + M 2 apless , (All) 

where M 2 apless is the (till unknown) contribution of the gapless fermions. A direct comparison with Eq. (|A5j) allows 
the identification 

f 2 = ^gapless ■ (A12) 

At this point, if we knew M 2 apless we should be able to determine f 2 in the low energy theory. However this 
contribution to the Debye mass is not well defined, as to calculate it one should choose an appropriate infrared cutoff 
and integrate out all quarks whose momenta are above such a cutoff. In this way, the theory would contain only the 
light degrees of freedom and the computation of M 2 apless would be straightforward. 

However, in the gapless regime we can use a different strategy: we already have computed f 2 using a completely 
different approach in the previous section, see Eq. (f2"5)) . Using this result, the matching condition (|A12[) can be used 
not to compute f 2 but to evaluate M 2 apless , namely 



M 2 gapless = ^-9(6^ - A) ^ A2 . (A13) 
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This result is equal to Mg R in Eq. (|A4[) . which is obtained by integrating in the loop over the quark momenta in the 
blocking region. This coincidence seems to show that the g2SC phase can be described like a two-component Fermi 
liquid, with a normal part given by the blocking region and relevant for the low energy dynamics, and a supcrfiuid 
part living in the pairing region. 

We consider now the Meissner mass: this case is instructive since it shows that the aforementioned interpretation 
of the g2SC as a two-component superfluid leads to meaningful results once one tries the matching. Following our 
interpretation we argue that the quarks in the blocking region, as they behave as a normal liquid, do not contribute 
to the Meissner mass. The same is true for the unpaired blue quarks. Finally, the quarks in the pairing regions are 
superfluid and do not contribute to the low energy dynamics. As a consequence, the only contribution to the Meissner 
mass in the low energy theory is given by the </> as in the 2SC phase, see Eqs. (|A8[) and (|A9|) . Matching with the high 
energy result in Eq. (|A6[) one gets 

^ 2 = ! # = m2 _ M2 m " _ M2 , (A14) 

J ' n D 1V1 gapless 1V1 Uue 

which leads to the result quoted in Eq. (|28|) once Eqs. (|A5j) . (|A6|) and (|A13j) are used. 
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